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Abstract. We study irreversible polymer adsorption from dilute solutions theoretically. Universal features 
of the resultant non-equilibrium layers are predicted. Two broad cases are considered, distinguished by the 
magnitude of the local monomer- surface sticking rate Q: chemisorption (very small Q) and physisorption 
(large Q). Early stages of layer formation entail single chain adsorption. While single chain physisorption 
times r a< js are typically micro to milli-seconds, for chemisorbing chains of N units we find experimentally 
accessible times r a d s = Q~ 1 N 3/ ' 5 , ranging from seconds to hours. We establish 3 chemisorption universality 
classes, determined by a critical contact exponent: zipping, accelerated zipping and homogeneous collapse. 
For dilute solutions, the mechanism is accelerated zipping: zipping propagates outwards from the first 
attachment, accelerated by occasional formation of large loops which nucleate further zipping. This leads 
to a transient distribution ui(s) ~ s~ 7,/5 of loop lengths s up to a maximum size s max « (Qt) 5//3 after time 
t. By times of order r a d s the entire chain is adsorbed. The outcome of the single chain adsorption episode 
is a monolayer of fully collapsed chains. Having only a few vacant sites to adsorb onto, late arriving chains 
form a diffuse outer layer. In a simple picture we find for both chemisorption and physisorption a final loop 
distribution f2(s) ~ s -11 ^ 5 and density profile c(z) ~ z~ 4 ^ 3 whose forms are the same as for equilibrium 
layers. In contrast to equilibrium layers, however, the statistical properties of a given chain depend on its 
adsorption time; the outer layer contains many classes of chain, each characterized by a different fraction 
of adsorbed monomers /. Consistent with strong physisorption experiments, we find the / values follow a 
distribution P(f) ~ /" 4/5 . 

PACS. 82.35.-x Polymers: properties; reactions; polymerization. 05. 40. -a Fluctuation phenomena, random 
processes, noise, and Brownian Motion. 68.08.-p Liquid-solid interfaces. 



1 Introduction 



Polymer layer formation is unavoidable when even weakly 
attractive surfaces come into contact with a polymer solu- 
tion |T|2*]. see fig.ffl Even for extremely dilute polymer so- 
lutions, polymer layers develop with densities which may 
be many orders of magnitude larger than the bulk poly- 
mer concentration 3 . This is due to the fact that by giv- 
ing up their bulk translational entropy, which costs a free 
energy of only kT, chains achieve an energy advantage 
proportional to the number of monomers per chain, N, 
for large enough N values. The topologically complex in- 
terfacial layers contain both surface-bound segments and 
large loops and tails extending into the bulk (see fig. [IJ. 
In principle, given enough time, adsorbed polymers are 
able to explore all accessible states ^01^] and the layer 
attains equilibrium. Equilibrium layers have been the fo- 
cus of a large body of experimental |3[71IH1IS|, theoreti- 
cal |2TU2ai23| work. 
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In a great many cases, however, the surface sticking en- 
ergy per monomer e exceeds fc^T. Available experimental 
evidence suggests that desorption processes and relaxation 
kinetics within the layer are then sharply slowed down 

Pll23ll3l23llIll23Gnl^ • 

For e values of several ksT time scales become so long 
that the layer build-up becomes essentially an irreversible 
process leading to non-equilibrium structure |.'S1II32| . One 
may speculate this is due to cooperative effects mediated 
by mutual topological chain constraints which drastically 
suppress mobilities near the surface jJUJ. Indeed, in the 
scaling theory 1 for equilibrium layers, when e exceeds 
fcgT the smallest loops in the layer become as small as 
the monomer size, a. Such small loops are likely to greatly 
restrict chain motion. 

For physical adsorption processes, large sticking ener- 
gies originate in hydrogen bonding, dispersion or dipolar 
forces or attractions between charged groups. Most metal 
and silicon-based surfaces are normally oxidized and many 
polymer species form strong hydrogen bonds with the sur- 
face oxygen or silanol groups. Biopolymers such as pro- 
teins and DNA attach tenaciously to many surfaces due 
to their large number of charged and polar groups |41I42| . 
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Fig. 1. The situation studied in this paper: polymer chains 
adsorbing onto a surface from a dilute polymer solution. We 
consider situations where monomer-surface bonds develop ir- 
reversibly, due either to chemical bonding or strong physical 
interactions, (a) Early stages of layer formation. The surface 
is almost empty and the first chains to arrive adsorb on the 
surface without interference from others, (b) At longer times a 
polymer layer of strongly interacting chains develops and chain 
densities on the surface become much higher than those in the 
bulk. 



Since hydrogen bonds, for instance, normally have ener- 
gies of 4fcsT and greater [E] it is apparent that strong 
physical bonds are very common. 

The extreme example of irreversibility arises in chemi- 
sorption 44 45 , 46 47114*%] where covalent surface-polymer 
bonds develop irreversibly (fig. [""}. In various technolo- 
gies polymers are attached by chemical reactions to solid 
surfaces either from solution as in colloid stabilization by 
chemically grafting polymers onto particle surfaces |49) . 
or from a polymer melt as occurs in the reinforcement of 
polymer-polymer [5(15 1152] or polymer-solid [53154155150] 
interfaces. In general, whether physical or chemical bond- 
ing is involved many applications prefer the strongest and 
most enduring interfaces possible and irreversible effects 
are probably the rule rather than the exception. 

This paper considers dilute solutions. Irreversible ad- 
sorption from melts and semidilute solutions, previously 
studied both experimentally |56M57llI^I59U6UM61l l"*"2 63 
and theoretically 64 65 66 , involves rather different pro- 
cesses. Our aim is to understand the structure and for- 
mation kinetics of layers which are formed under these 
irreversible circumstances where the usual statistical me- 



chanical approach is inapplicable to the non-equilibrium 
structures which form. A series of experiments by Granick 
and coworkers [23 E3 EHI EH E2I EDI E31 EH ESI have ex- 
amined these issues. These workers found that when e 
reached values of only a few kT, polymer relaxation times 
became large and equilibrium layers were not attained. 
This was most clearly apparent in experiments follow- 
ing polymethylmethacrylate (PMMA) adsorption onto ox- 
idized silicon (e 4fceT) via hydrogen bonding from a di- 
lute CC1 4 solution [30133131132] . Measuring both the total 
adsorbed mass, r, and the surface-bound part, /bound, as 
a function of time by infrared absorption spectra, it was 
found that early arriving chains had much higher fractions 
of bound monomers, /, than late arrivers, and these / val- 
ues were frozen in throughout the experiment's duration of 
several hours. Essentially monomers remained irreversibly 
fixed to the site they originally landed on, or at least close 
to this site. Measuring the distribution of / values among 
chains they found a bimodal distribution shown in fig.[2Ib) 
with two peaks at small and large /, respectively. This is 
strikingly different to equilibrium layers where all chains 
within the layer are statistically identical and are charac- 
terized by the same time-averaged value, / = ^bound/^°°' 
where oo denotes asymptotic values (t — > oo). 

A number of analytical and numerical efforts [4611471 
|<37J|fDHl"SSllZ01l7^IS2 have addressed irreversible polymer 
adsorption from solution. However, the understanding of 
these phenomena remains very far from the quantitative 
level which has been achieved for equilibrium layers. In 
this paper we develop a theory which amongst other objec- 
tives aims to understand the experimental findings of refs. 
[30133131132] . We emphasize adsorption from a dilute poly- 
mer solution under good solvent conditions and consider 
systems where monomer-surface bonds are strong enough 
that they are effectively irreversible within the experimen- 
tal timescales. Certain results are generalized to theta sol- 
vent solutions. We calculate the relationship IboundC-O, 
how each of /bound and r depends on time, as well as the 
final layer's distribution of chain contact fractions, P(f). 
We will compare our predicted final layer structure to the 
well established theoretical results for equilibrium layers 
which predict a density profile c(z) ~ z~ 4 / 3 , and a self- 
similar loop distribution, i?(s) ~ s~ n / 5 . Similar to the 
picture that was developed by Granick and co-workers, we 
find that the final layer consists of two populations: early 
arriving chains lie flat on the surface while late arrivers 
can only adsorb with an ever-decreasing fraction of their 
monomers onto the surface leading to a diffuse weakly 
attached outer layer. We find a universal distribution of 
/ values, P(f) ~ / -4 / 5 for small / which agrees rather 
closely with the experiments of refs. [31132] . Interestingly, 
we find layer loop distributions and density profiles obey- 
ing the same scaling laws as those of equilibrium layers. 
Chain configurations are very different, however, leading 
to radically different physical properties of the layer. 

Our picture of irreversible polymer adsorption is in 
some respects qualitatively similar to the theoretical one 
of the workers of ref. who simulated their experiments 
in a random sequential adsorption [72j framework. They 
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Fig. 2. (a) Theoretically predicted dependence of surface- 
bound mass, /bound, on total adsorbed mass per unit area, F. 
For irreversible physisorption an initial linear regime crosses 
over to sharp saturation Zi/b oun d ~ (Ar) 6 as asymptotic cov- 
erages are approached. A similar curve has been measured ex- 
perimentally in ref. 13 1 1. For chemisorption the relation is the 
same, but with an additional initial regime /bound ~ r s ^ 3 . (b) 
Probability distribution of fraction of bound monomers / per 
chain. Grey bars reproduce experimental data from ref. [31] . 
Empty bars are the theoretically predicted P(f) ~ / -4 ^ 5 for 
/min < / < u where values for / m i n = 0.09 and ui — 0.47 were 
derived from the measurements in ref. |31| . The delta function 
at / = lo is the contribution from early arriving chains and is 
expected to be broadened in reality. 



visualized chains as "deformable droplets" : at the early 
stages when the coils arrive onto a bare surface, each 
droplet adsorbs a certain maximum cross-sectional area. 
As available surface area for adsorption become scarce, in 
order for late-arriving chains to adsorb, it was assumed 
droplets deform by reducing their cross-sectional areas 
parallel to the surface to fit into the empty space. In so 
doing, they become more extended into the bulk. Using 
this model they generated a P(f) similar to the experi- 
mental one of fig. |2Jb) ■ The picture developed here differs 
in this respect: late arriving chains freely overhang early 
flat-lying chains and rather than fit into a single available 
connected surface area they attach at disconnected empty 
surface sites. 




Fig. 3. The two classes studied in this paper: physisorption 
and chemisorption. Simplified view of monomer free energy as 
a function of distance between monomer and surface. In ph- 
ysisorption there is no activation barrier and monomer-surface 
association is immediate upon contact. When sticking energy 
e exceeds a few kT, experiment indicates that effective des- 
orption rates become very small, presumably due to complex 
many-chain effects. Chemisorption typically involves a large ac- 
tivation barrier, u ^> kT, which needs many monomer-surface 
collisions to traverse. The adsorption strength is also large, 
e S> kT. Some systems are in practice mixtures of chemi and 
physisorption, a complexity not dealt with in the present work. 



The process of irreversible polymer layer formation 
entails progressive freezing in of constraints due to irre- 
versible monomer-surface bonding. These constraints grad- 
ually reduce the volume of configurational phase space 
available. Thus ergodicity is inapplicable and in lieu of 
the algorithms of equilibrium statistical mechanics one 
must follow the kinetics of chain adsorption and layer 
build up. It is important to distinguish carefully between 
two broad classes of adsorption kinetics, physisorption and 
chemisorption, which are characterized by very different 
values of the local reaction rate Q. We define this as 
the conditional monomer-surface reaction rate, given an 
unattached monomer contacts the surface (see fig. |!IJ . In 
physisorption, monomer attachment is essentially diffusion- 
limited, Q « l/t a , where t a is monomer relaxation time. 
Chemisorption, where adsorbed monomers form very strong 
chemical bonds with the surface |Y3. 44. 45. 14611471148] is 
much slower 74 ,52 ,75 with Q values 8 or more orders of 
magnitude smaller than those of physisorption. The origin 
of this difference is that chemical bond formation usually 
involves a large activation barrier (fig. OH . 

Our theoretical approach is to make the simple as- 
sumption of total irreversibility, motivated by the exper- 
iment discussed above: once a monomer bonds with the 
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surface, our model is that this bond never breaks. The 
processes of chemisorption and physisorption will be an- 
alyzed separately. Though we find both lead to similar 
final layer structures, the kinetics are very different. In 
physisorption a single chain adsorbs onto a bare surface 
in a time of order the coil relaxation time or less, typ- 
ically microseconds, whereas single chain chemisorption 
may last minutes or even hours and is thus observable 
experimentally. We begin with a study of chemisorption 
in sections 2, 3, and 4. Specifically, in Section 2 we show 
that quite generally there are three possible modes of sin- 
gle chain adsorption, depending on the value of a certain 
critical exponent 9. Single chain chemisorption from good 
and theta solvents is then studied in detail in Section 3. In 
Section 4 we consider the later stages of the kinetics when 
chains overlap and dense layers are formed (fig. . The 
case of irreversible physisorption is treated in Section 5. 
We compare our results to experiment and conclude with 
a discussion in Section 6. An announcement of our results 
has appeared in ref. |76| . 



2 Chemisorption, General Discussion: Three 
Modes of Adsorption 

Our interest is an initially bare surface confronting a di- 
lute solution of functionalized chains. During the earlier 
stages of surface layer formation, the coverage is small and 
individual chains do not see each other. In this section we 
consider in detail how a single chain chemisorbs onto an 
interface, remembering that chemisorption is character- 
ized by very small values of the local monomer-surface 
reaction rate Q. The chain will make an initial attach- 
ment and then develop a certain loop profile as successive 
monomers gradually attach, eventually leading to com- 
plete chain adsorption in a certain time r a( j s (see fig. 
whose dependence on chain length is an important char- 
acteristic. Since the early polymer layer consists of chains 
dilute on the surface as in fig.^a), the initial layer struc- 
ture is a superposition of these single chain loop profiles. 

In this section we consider chemisorption in its full 
generality. We show that there are three distinct classes 
of behavior, each with different adsorption modes, loop 
profiles and adsorption times r a d s . Which class a given 
experimental system belongs to depends on the bulk con- 
centration regime and other factors. The classes are de- 
fined by the value of a critical exponent characterizing 
polymer statistics near surfaces. 

We will assume that all monomers are identical and 
chemically active, and that the surface has a free ener- 
getic preference for the solvent over the polymer species. 
This means that in terms of physical interactions the poly- 
mers see a hard, repulsive wall. We choose units where the 
monomer size a is unity. 

Consider the chain in fig. 0Jb) whose first monomer 
has just reacted with the surface. Which of the chain's re- 
maining monomers will react next? What is the sequence 
of chemisorption events? This depends on the form of the 
reaction rate k(s) at which a reactive group s monomer 



units along the backbone away from the bound monomer 
reacts with the surface, as shown in fig. Ufa). Now the im- 
portant feature of chemisorption is that due to the small- 
ness of Q in order for any chemical bond to form a time 
much larger than the attached chain's relaxation time is 
needed. Thus quite generally the reaction rate for a given 
monomer at any moment is proportional to the equilib- 
rium contact probability of finding this monomer on the 
surface given the current constraints due to all chemical 
bonds formed at earlier times. In this particular situation, 
this is 

*i«) = ff, a) 

where Z{N) is the partition function of the chain with one 
monomer bound (middle of fig. Ufa)) and Z(s, N) is the 
partition function of the chain which has the additional 
constraint that the s th monomer is also bound (last of 
fig. Ufa)). Physically one expects that for small enough 
s, this is independent of N, i.e. p(s\N) « s~ e where the 
value of the exponent 9 reflects the equilibrium polymer 
statistics, given the constraints. This then leads to the 
following expression for the reaction rate: 



k(s) 



Q 



(a<N) 



(2) 



The total reaction rate for the next adsorption event is a 
sum over all the N — 1 monomers which may adsorb next. 
These belong to the two tails in fig. Bib) which are of order 
N in length. The net reaction rate is thus approximately 



JY 



k(s)ds 



(3) 



which, depending on the value of 9, is dominated by either 
small or large s. 

This argument is then repeated for the reactions of 
the second, third and subsequent monomers, in all cases 
described by a rate with the same small-s behavior as in 
eq. (O). The only difference is that tails are replaced by 
loops. Thus the entire kinetic sequence is characterized by 
the single exponent 9. The nature of the kinetics depends 
on the value of 9 as follows. 

(i) 9 > 2. Suppose the equilibrium statistics (eq. QJ) 
are such that k(s) decays faster than 1/s 2 . Then T^totai ~ 
k(l) and typically a monomer near the first attachment is 
most likely to attach next. The third monomer to attach, 
repeating the same argument, will be near the first two, 
and so on. In this case therefore the chain would zip onto 
the surface starting at the first attachment point, as shown 
in the top of fig. Etb) . Since each new attachment occurs 
at the same rate each time, the full zipping time would 
then be r ads w N/k(l) Q^N. 

(ii) 9 < 1. In this case T^totai is dominated by s of 
order N reflecting the fact that even though a monomer 
with small s is on average more frequently near the sur- 
face than a distant one, there are many more distant ones 
and their number is the dominating factor. Thus a distant 
monomer of order N units away from the first graft at- 
taches next, leading to the formation of a loop of size of 
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Fig. 4. Schematic of single chain chemisorption. (a) A chain 
makes its first attachment with the surface and starts its ad- 
sorption process. The rate of first monomer attachment is pro- 
portional to Q, namely the reaction rate given the monomer 
touches the surface. Subsequent monomer adsorption is deter- 
mined by the form of the reaction rate k(s) for the s th monomer 
away from the first attachment. For s < iV, where N is length 
of the tail, k ~ Q/s e where the exponent 9 reflects equilib- 
rium polymer statistics near interfaces, given the constraints 
imposed by earlier reactions, (b) Depending on 9, three modes 
of adsorption are possible, (i) Zipping, 9 > 2. Monomers ad- 
sorb sequentially forming a long train, (ii) Accelerated zipping, 
I < 9 < 2. As (i) but zipping is accelerated by nucleation of 
new zipping centers due to formation of large loops, (iii) Uni- 
form collapse (9 < 1). Monomers distant from adsorbed ones 
come down first leading to formation of large loops whose size 
is uniformly decreasing with time. Chemisorption from dilute 
solution with good solvent is accelerated zipping (9 = 8/5). 



zipping since large loop formation effectively short-circuits 
the pure zipping process by nucleating new sources of zip- 
ping as shown in the middle case of fig- Bib). Now eq. (QJ 
implies (unlike 9 < 1) larger loops take longer to form. 
Thus by tn loops of all sizes have come down, i. e. the 
chain has adsorbed and we conclude r a d s ~ tn. 

In summary, depending on the value of 9 the chemisorp- 
tion time has three possible forms: 

{const. (9 < 1) 
N 9 ' 1 (1<0<2) (5) 
N (9>2) 

In the next section we calculate the value of 9 under 
good solvent conditions and find that it belongs to the 
accelerated zipping case. We then study the corresponding 
case in detail. 



3 Single Chain Chemisorption in a Good 
Solvent: Accelerated Zipping 

In this section we consider the kinetics of single-chain 
chemisorption, describing polymer layers during the early 
stages of the chemisorption process, fig.^a). Expressions 
are derived for the number of bound monomers, 7bound(T), 
and the single-chain loop distribution profile, cj r (s), where 
r measures time after first monomer attachment onto the 
surface. In order to perform these calculations, in this 
section we first evaluate the reaction exponent 9 of eqs. 
(|2j and and find it belongs to the accelerated zipping 
class. Subsequently we solve the accelerated zipping kinet- 
ics, first using simple scaling arguments and then a more 
detailed solution of the rate equations. 



order N. Repeating the argument for the subsequent re- 
actions leads to a process of uniform collapse of the chain 
where roughly centrally positioned monomers of loops re- 
act with the surface at each step, as shown in the bottom 
case of fig. 0] In this case the average lifetime of a loop 
with s monomers is Ti 00 p(s) ~ 1/ ds'k(s') ks Q~ 1 s 1 ~ 9 . 
Thus smaller loops take longer time to collapse and the 
rate limiting step for full adsorption is the collapse of the 
smallest loops: r ads w r loop (l) w Q' 1 . 

(iii) 1 < 9 < 2. This is intermediate between zipping 
and collapse. Though T^total is dominated by small s, sug- 
gesting zipping, there is in fact enough time before pure 
zipping completes for large loops to form. To see this, re- 
turn to the 1-graft situation on the extreme left of fig-GI D ). 
The time for a loop of order N (for argument's sake, big- 
ger than N/2) to form from one of the two tails extending 
from the grafted monomer is roughly: 




N/2 



ds k(s) 



(4) 



This is much smaller than r a( j s ~ N which pure zipping 
would give, and thus before puring zipping is complete, 
large loops must have formed. We call this case accelerated 



3.1 The Reaction Exponent 9 

To our knowledge 9 in good solvents has not been cal- 
culated before. We show here that it can be expressed 
in terms of other known critical exponents characterizing 
polymer networks |77U78| . The simplest way to derive 9 is 
to consider a loop of N monomers bound to the surface 
by its two ends as in fig. |S[a). Then the reaction rate of 
its s th monomer resulting in the formation of two loops of 
lengths s and N — s, as shown in fig. EJb), is 



k(s\N) = Q 



Zioopjs, N - s) 
Zio OP (N) 



(6) 



where Z\ oop (N) and Z\ oop (s, N — s) are the partition func- 
tions of the loops of fig. EJa) and[SJb), respectively. Now 
although Zi oop (N) is known, Z\ oov (s, N) is only known for 
s of order N. From ref. 1771: 



7 loop 



(N/2, N/2) w n N % 



Z\oop(N) w [± N la ~ 

where fi is a constant of order unity and the values of of 
7o,7c have been calculated analytically and numerically 
[77|ES1 • Now demanding that Zi oop (s,N — s) follows a 



6 



Ben O'Shaughnessy, Dimitrios Vavylonis: Irreversible Adsorption from Dilute Polymer Solutions 




Fig. 5. A polymer loop anchored to the surface by its two 
ends. The reaction rate of its s th monomer with the surface, 
leading to two loops of length s and N ~ s, is proportional to 
the ratio of the partition functions of the loop before and after 
the reaction. 



power law for small s and that it reduces to eqs. (JZJ) in the 
s — > 1 and s — > N/2 limits one has 



Z loop (s,N- S )*n N N^-\(jj 



Here the scaling function £ must obey 



x-° (x -> 0) 

1 (a; =1/2) 

{l-x)- e (x -> 1) 



(8) 



(9) 



where the x — > and x — > 1 behaviors are identical by 
symmetry and 



# = l + ^«8/5 (good solvent) 



(10) 



which follows after using the identity [7? j a — 7c = v + 1. 
Here i/ « 3/5 is the Flory exponent [7J3. From eqs. ©, 
0, and |JS} one thus has 



(11) 



which for s <C iV, reduces to eq. (J2J, fc ~ s -9 . Thus single 
chain chemisorption in a good solvent is characterized by 
8 k 8/5; since 1 < 8 < 2 the adsorption mode is acceler- 
ated zipping. 

Eq. was derived for the particular case of a loop as 
in fig. El But even if the reacting monomer is part of a tail 
as in fig. P| a), by a straightforward generalization of the 
reasoning of this section using the exponents calculated in 
ref. [77], one can show that the reaction rate has the same 
form for s <C N. One can show the same is true if the ends 
of the loop in fig. are connected to other loops. 

The above results for 8 directly generalize to theta sol- 
vent solutions. In this case polymer statistics are effec- 
tively ideal [T^j and v — > 1/2 in all of the above, i. e. 



8 = 3/2 



(theta solvent) 



(12) 



which also corresponds to accelerated zipping. In fact, 
for this case the exponent can be obtained more simply as 
follows: the probability that the s th monomer, measured 
from a given surface attachment, is in contact with a hard 
wall is proportional to the return probability P of a ran- 
dom walk which starts one step away from an absorbing 
surface after s steps, i. e. k ~ P ~ s 



Finally we remark that the partition functions of eq. 
correspond to loops such as those of fig. whose end 
locations on the surface are annealed [77]. However for 
chemisorption onto solid surfaces the ends are either com- 
pletely fixed or may diffuse very slowly on the surface. 
(Note that chemisorption onto liquid interfaces are cases 
where the end locations are truly annealed.) Nonetheless, 
all our results for s <C N and the scaling structure of eq. 
(|llfl must remain the same even in this case since in the 
s ^ N limit, k(s\N) is independent of the location of the 
other end. Our general conclusions are thus expected to 
be robust, provided the two ends are not so far apart that 
there are strong lateral loop stretching effects. In the fol- 
lowing we self-consistently assume that as the adsorption 
process proceeds there is no tendency for generated loops 
to be in such stretched configurations. 



3.2 Accelerated Zipping: Scaling Analysis 

Let us consider now the kinetics of adsorption starting 
from a polymer chain as in fig.QIa) which has just made its 
first attachment with the surface with an interior monomer. 
This is in fact typical: we show in Appendix A that chains 
are much more likely (by a factor N 0A8 ) to make their first 
surface contact with an interior monomer because there 
are many more (of order N) such monomers as compared 
to chain ends. 

Now the chain starts to chemisorb at r = 0. From the 
previous subsection we have seen that in the accelerated 
zipping mode the chain adsorption time T a( js is equal to tjv , 
the timescale associated with loops of size N. What will 
the chain's loop profile be for times r smaller than r^r? Ini- 
tially, since the reaction rate is dominated by small s, the 
chain will start to zip and sequences of bound monomers 
("trains") will grow outwards from the first attachment 
point. As time proceeds though, large loops start to come 
down due to adsorption of monomers distant to the al- 
ready bound ones. These loops have a certain distribution 
of loop sizes ui T (s) and their effect is to accelerate the 
zipping process since they nucleate further sources of zip- 
ping and train formation. The maximum size of such loops 
which had enough time to form by time r, i. e. for which 



f s J L x ds'k(s') a 1/r, is 



(Q- 



(13) 



-3/2 



which increases with time. The maximum loop size thus 
becomes of order N sXtm- 

Overall, we expect the chain to consist of three parts, 
shown in fig. HJ] (i) Two large tails, each of length of order 
N . These are the two tails of the polymer chain which was 
initially bound by one of its interior monomers as in fig. 
Eta). For t <^Tjy the tails had neither enough time to de- 
cay into large loops, nor to completely zip. (ii) Loops with 
a loop distribution lo t (s). (iii) Trains of bound monomers 
whose number we define to be 7bound(7")- 

Now 7bound('?") is easily found by making the ansatz 
that it follows a power law in time. In addition, for short 
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Fig. 6. Snapshot of loop profile of a chemisorbing chain dur- 
ing its accelerated zipping down onto a surface. At an earlier 
time the chain adsorbed with an interior monomer from which 
zipping started to propagate outwards accompanied by occa- 
sional formation of large loops which nucleated further zipping 
sources as shown. The distribution of loop sizes follows a uni- 
versal power law and as larger loops form with increasing time, 
loop sizes are larger near the outermost tails. The tails act as 
a continuous loop source thus sustaining the loop distribution. 
Eventually, after r a d s ~ Q~ 1 N V the largest loop becomes of 
order the tail size, N, the tails are completely consumed and 
the structure collapses onto a completely flat configuration. 



times 7bound is independent of N; one can imagine send- 
ing the chain size to infinity, which would not affect the 
accelerated zipping propagating outwards from the initial 
graft point. This dictates: 



7bound(T) w N 



TN 



1/u 



(r < t n ) (14) 



since by tn most of the chain has adsorbed, i. e. we must 
have 7bound(Tjv) ~ N. 

We now evaluate u T (s) by making the ansatz that it 
also has power law behavior: 



LU T (s) 



1 



(1 « a <C s" ax ) (15) 



with 6 > 1. Since by r there has been enough time for 



the formation of order one loop of length 



the nor- 



malization has been chosen such that there is one loop 
of order s™ ax in size. Now the total number of loops is 
thus L(t) = j x dsuj T (s) ps (s" lax ) 5_1 . These loops pro- 
vide L(t) new nucleation points for further zipping. The 
reaction rate at each such point is dominated by small s 
(since 9 > 1) and is thus of order fc(l), so 



^7b 



dr 



k(l)L(r) 



,(S-l)/u 



(t<tat) . (16) 



Now demanding that (|16fl reproduces eq. (|14[> the value 
of S is determined: 



5 = 2-1/ 



(17) 



Notice that the above results are self-consistent: the num- 
ber of loops L(t) is much smaller than 7bound('?"), and 



essentially all reacted monomers do belong to trains as 
assumed in eq. (|14|) . 

Let us summarize this subsection's results, setting v — 
3/5 which is the value relevant to our main interest of 
dilute solutions with good solvent. From eqs. I|13|) . (|14|l 
and H15|) we have: 



7bound(T) W (Qt) 



5/3 



valid for r < rjy and s < s" 1. 



lu t (s) ps (Qt) 



2/3-7/5 



(18) 



3.3 Accelerated Zipping: Solution of Rate Equations 

In this section we analyze the accelerated zipping phe- 
nomenon in a more rigorous way by solving the loop evo- 
lution dynamics. We will recover all of the scaling results 
of section 3.2. The time evolution of the chain's loop dis- 
tribution is given by |46j 



duj T (s) 
dr 



N 



= 2 J ds' uv(s') k(s\s') - J ds'uj T (s)k(s'\s) , 

where the first term on the right hand side (rhs) describes 
the rate of formation of loops of length s by bigger ones, 
while the second term on the rhs is the rate of decay of a 
loop of length s into smaller loops SO and k is given by 
eq. (flip . The initial condition for eq. I|19[l is wo(s) = 6(s — 
N), i. e. there is only one initial loop of length N (in the 
following for simplicity we do not distinguish between tails 
and loops). A basic assumption in eq. Ijl9|l is that k(s\N) 
retains the form of eq. throughout the adsorption 

process. That is, we assume topological and many-loop 
excluded volume effects (beyond those contained in eq. 
©) are not strong enough to alter the essential form of 
k(s\N). It is easy to show by integrating eq. (|T§|) over 
all s (see calculation in Appendix B) that the kinetics 
conserve the total number of monomers. Formally, uj t (s) in 
eq. (|19J) is an ensemble average over many chains. However 
the number of loops formed by a single chain becomes 
much larger than unity as time increases and thus eq. (|19|l 
accurately describes the loop distribution of a single chain 
as well. 

Although both integrals in eq. I|19f) diverge at s' = s 
and s' = 0, these divergences cancel with one another. In 
fact defining D T (s) to be the number of loops larger than 



rfs'w r (s') 



(20) 



we obtain the following equivalent rate equation which 
has well-defined integrals: 



dD T {s) 
dr 



N 



ds' \{s\s')u T {s') , 



(21) 



where 



o r 1 - 3 / 3 ' 

\(s\a') = ±L dxC(x) 

s Js/s' 



1-s/s' 
s I s' 

/ -AQ/(s' 
\BQ/s" 



{af - a) 
(s' » a) 



(22) 
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Fig. 7. Sketch of the size probability distribution ui T (s) for 
the large tail, (s — * N, figEJ for a polymer chemisorbing from a 
good solvent, see eq. 1241 . The initial (r = 0) tail length is a S- 
function peaked at N and decreases with time r on average by 
s™. The broad distribution at r represents fluctuations in tail 
length among a population of many chains. The normalization 
of the distribution remains of order unity for r -C tn- 



with A, B positive constants of order unity. Eq. Q21|l 
is derived by integration of eq. I|19l) . Its validity is easy 
to check by differentiation with respect to s and using 
Q{s/s') = C(l — s /s')] one then recovers eq. lfT§| . Notice 
that the function A(s|s'), which describes how much each 
loop s' contributes to changes in D T (s), is negative for 
s' < 2s since a loop shorter than 2s always creates at 
least one loop smaller than s. Also A(s|s') is positive for 
s' > 2s since at least one loop bigger than s is created by 
a loop longer than 2s. 

We saw in subsection 3.2 that the two large tails form- 
ing after the first attachment have of order N monomers 
until essentially the end of the chemisorption process. In 
order to describe the tail kinetics, let us consider s very 
close to N in eq. l(2l|l: 



dD T (s) 
dr 



N AO 

ds' , . uj t (s') (s->N) . (23) 



(s' - s) v 



Self-consistently, one can show that error terms arising 
from approximating A by —AQ/(s' — s) v in eq. i|23[) are 
higher order. Eq. I|23fl is solved in Appendix C, with solu- 
tion 



uj t (s) = 



where 



1 



..max \ _/V — S 

h{x) - 



l+v 



N - s 



(x > 1) 
const, (i « 1) 



(N - s <C N) 
(24) 



The cutoff function h(x) decays to zero exponentially 
rapidly as x — > oo. Eq. l(2^|) (see fig. [7J describes the 
length of the two initial tails, which is continuously de- 
creasing. After time r the length has decreased by s™ ax . 
This decrease is the maximum loop size that any tail could 
have created by time r. Due to fluctuations, the initial 5- 
function peak broadens as it moves towards s = 0. Thus 
w r (s) in eq. (|24|l is the probability distribution of the tail 



CO 




Fig. 8. Log-Log plot of single-chain loop distribution lu t (s) for 
chemisorption from dilute solutions at two different times, ri < 
T2- The distribution consists of an inner and an outer region 
with to ~ s~ 7,/5 and u ~ s _8//5 , respectively, for good solvents 
(y — 3/5). Both the normalization which is equal to the total 
number of loops, and the size of the inner region increase with 



time. The outer region extends up to s of order N 



For 



larger s the distribution is shown in fig. |7| describing the size 
of the large tail. 

length. One sees that by time Tjy the length of the tail (the 
largest non-adsorbed chain segment) shrinks to zero, thus 
marking the end of the single chain adsorption process. 

Consider now loops (1 <C s <C N) generated by oc- 
casional grafting of distant monomers. For short enough 
times the only sources of loop formation are the tails which 
dominate the integral on the rhs of eq. I|21(l at s' of order 
N. Assuming for the moment that this continues to be 
true for all times, one has from eq. (|21|) 



dD T (s) Q 
dr ~ s v ' 
which after integration gives 



LO T (s) 



1 



l+u 



( S ; nax « s « N) (25) 



(sr«s«^) • (26) 



Thus the first loops which form have the same loop distri- 
bution profile as the small s behavior of k(s\N). However, 
the validity of eq. is limited to s > s™ ax . This can 
be seen by substituting back eq. H26(l into eq. H21(l . reveal- 
ing that the assumption that only tails contribute to the 
rhs of eq. I|21l) is incorrect for s ~ s" lax . This limits the 
validity of eq. f2^|l to those loop sizes too great to have 
been formed by tail collapse events by the time r. This 
is reflected by the fact that the integral of eq. (2G) in its 
region of validity is of order unity and dominated by s of 
order s™ ax . 

Now for s < s™ ax , the main body of the loop distri- 
bution, we seek a quasi-static solution, dD T /dr ~ 0, or 
equivalently from eq. (|21|) . 



N 



ds' X(s\s')uj t {s') w 



(1 « s < s™ 



(27) 



representing the expectation that such loops reach a self- 
maintained universal distribution. We show in Appendix 
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D that eq. 127J1 has a power law solution: 



w T (s) 



s™ ax ) . (28) 



The normalization of eq. ((28(1 was fixed by demanding 
continuity with eq. 1)26(1 at s" lax . Combining eq. 1)28(1 and 
eq. ((26(1 we obtain the overall loop distribution shown in 

fig-El 

Now the remaining part of the loop profile is the partic- 
ular contribution due to trains. That is, we must still cal- 
culate the total number of adsorbed chain units 7bound(T)- 
One might think this is just w(l), the number of minimal 
length loops. This is wrong because there is in fact a sink 
for loops of length 1, i. c. a current out of the uj(s) distri- 
bution into the total mass of trains, 7bound('?")- The latter 
obeys the dynamics of eq. 1(19(1 but with the decay term 
deleted since bound monomers cannot decay into smaller 
ones: 



ound 



dr 



ds'uj T (s')k(s\s r ) sa Q" 1 L(r) . (29) 



Here unlike eq. ((19(1 we used an explicit cutoff at s = 1. 
In the last expression of eq. lj2"9l we used the definition 
of the total number of loops, L(r), and took into account 
that the integral is dominated by its lower limit. Thus eq. 
l{2"§ll reduces to eq. (tTB)l leading to the solution of eq. (Till) . 

In summary, using v — 3/5, our results for single chain 
chemisorption are 



7bound(T") 



(Qt) 5 / 3 
N 



(T < T N ) 
(T > T N ) 



as in fig. ^ Initially, this is a simple superposition of sin- 
gle chain adsorption processes discussed in section 3. As 
chains build up and overlap, this becomes a many-chain 
phenomenon. We will determine the kinetics of total and 
surface-bound mass per unit area, r(t), and /bound (i), 
the relationship /bound (r), the monomer density profile 
c(z), and the internal structure of the layer characterized 
by the loop distribution per site, j?(s). Finally we deter- 
mine the distribution of fraction of adsorbed monomers 
per chain, P(f). We explicitly set v = 3/5 and it will be 
convenient to reinstate explicit reference to the monomer 
size a, hitherto set to unity. 



4.1 Monolayer Forms at Early Stages 

In the early stages of adsorption, surface-attached chains 
are dilute on the surface and each one performs its ac- 
celerated zipping down in a time tjv- Thus the layer is a 
superposition of such chains which arrived on the surface 
at different times, the structure of each being given by 
eq. (130(1 . Now the rate of chain arrival onto the surface 
is a 2 dr/dt = QNtfi sur {, where _T is monomers per unit 
area and </> su rf is the equilibrium monomer density at the 
surface: 



^surf 



fbulk" 



Z SU rl(N) _ 4>- 



Jbulk 



(iV) 



bulk 

aT 



(32) 



Here Z SUI { , ^buik are the single chain partition functions 
given a monomer on the surface or in the bulk, respec- 
tively. Their ratio has been shown to be 1/N in ref. [77]. 
Solving for the kinetics one thus obtains 



(s™ a 
(s™ a 



r{t)a 2 = ^buikQt 



t < t: 



chem 



l/(Q<j> 



bulk) 



(33) 



x /s) 7/5 
sc/ s )8/5 

7 (AT - s)] 8 / 5 h(sf ax /(N - s)) 



(1< s < s™ ax ) 
(s™ 3 * <s<N) 
(s^N) 

(30) 



where 



r N « Q _1 A 3/5 



(Q T ) 



5/3 



(31) 



and the expression for uj t is for t < t^. A schematic of 
the loop distribution is shown in figs. El and |H1 Apart 
from extra details these results are identical to the ones 
we obtained in section 3.2, eq. l(TB|) . using scaling argu- 
ments. Analogous results to eq. ((30(1 can be obtained for 
theta solvents replacing v — > 1/2 in all expressions of this 
section. 



4 Chemisorption: Kinetics of Polymer Layer 
Formation 

In this section we consider the full chemisorption process 
where many chains simultaneously attach to the surface 



where t^ m is the timescale at which -T « 1/a 2 , i- e. when 
surface density is of order one monomer per "site" of area 
a 2 . Clearly, as t^ m is approached, the kinetics of eq. (|3*3*|l 
should be drastically modified to take into account the 
depletion of available landing sites for new chains. 

Notice that r(r N )/N w (fauik/ 4>*) / R% where <t>* = 
jy-4/5 j g cnam overlap threshold concentration |79j. 
Thus for dilute solutions, </>buik < 0*, attached chains are 
dilute on the surface at time tat, the necessary time for 
a chain to fully adsorb. At any moment, therefore, those 
chains actually in the process of chemisorbing are dilute 
on the surface and so do not interfere with one anothers' 
adsorption process. Chains adsorb independently. (Note 
that for higher concentrations, i. e. for 0buik in the semidi- 
lute regime, chain interference effects would be important 
and the resulting layers different.) 

Although a given adsorbing chain does not see simulta- 
neously adosrbing chains, in the later stages of layer forma- 
tion it will be affected by those which have previously ad- 
sorbed, because these diminish the available empty surface 
sites. This effect becomes strong after time ^at™- Some in- 
terference of this type will commence at an earlier stage, 
when flattened down adsorbed chains first start to overlap 
on the surface. This happens at coverage r — N/R 2 F « 
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N~ x f 5 a 2 , which occurs after time i ove ria P ~ t^W" 1 / 5 . 
This is initially a very small effect since at time Overlap the 
fraction of unavailable sites is small for large N, of order 
TV -1 / 5 , so the essential features of the adsorption kinet- 
ics will remain unmodified. This becomes an order unity 
effect only at t^ m . Between the two times there is a con- 
tinuously increasing interference. However, for practical 
values of N this is such a weak power that the timescales 
^overlap and tg~g m will be difficult to distinguish. (In fact 
for theta-solvents, v — > 1/2, the two timescales become 
the same). 

Now the bound component of the attached chains of 
eq. (|33|1 is a sum over the bound mass of each chain which 
depends on the time it arrived on the surface: /bound (t) = 
J*dtT(t')j houlld (t - t')/N. Using eqs. (33J and (jUJ one 
thus has 



Abound (t) a 2 ~ 

^buikQiV 3 / 5 ^/^) 8 / 3 (t < T N ) 

4>bulkQt (t n < t 



^ j-chcni\ 
^ "-sat I 



(34) 



The short-time kinetics of eqs. (|33f) and l|34(l are sche- 
matically shown in fig. 1111 One sees two regimes in the 
short time behavior of /bound which are also reflected in 
two regimes of /bound (r). From eqs. l|3^|) and (|3"4*|) one has 

r(r N )[r/r(r N )} 8 / 3 , r<r( TN ) 
u>r, r( TN ) <r < a - 2 

This early portion of the /bound (-0 curve is illustrated 
in fig. [Ha). In eq. I|35(l we explicitly introduced the pro- 
portionality prefactor u>. This is a small-scale species- 
dependent constant representing the fact that even for iso- 
lated chains, steric constraints at the monomer level pre- 
vent every mono-mer from bounding. Thus u> is the frac- 
tion of monomers which are allowed by such constraints 
to bound. 

In summary, during the early stages of chemisorption, 
chains flatten out on the surface uninhibited by the pres- 
ence of others. For very short times, t < rjy, none of 
the attached chains has completed its adsorption and the 
layer's loop distribution is a superposition of single-chain 
loop structures given by eq. (j30Jl ■ summed over different 
arrival times t — r. For times longer than r^r, essentially 
all attached chains have fully adsorbed and a monolayer of 
flattened chains starts to develop, which at times of order 
tf^ m has almost covered the surface. 

Up to igat m the fraction / of bound monomers is ap- 
proximately the same for all chains and equal to to. Thus 
P(f) during this stage is sharply peaked at to. In reality, 
we expect two types of effects will somewhat broaden this 
distribution. The first is fluctuations in / values around 
u due to random events typical of multiplicative random 
processes characterizing irreversibility (e. g. monomers tra- 
pped in knots might not be able to bound to the surface). 
Such fluctuations would be interesting to characterize nu- 
merically following the example of ref. in which a 
monte-carlo method was used to study the structure of 
a single fully collapsed chain on a surface and an u> value 




Fig. 9. Typical configuration of adsorbed chain at late stages 
of adsorption. Such chains can adsorb only onto free empty 
sites (shown as gray discs) which are separated by l scp . They 
thus form loops of s monomers, with as 3//5 = l scp . 



was extracted using 7bound ~ ojN for large N. We antici- 
pate a second source of broadening for chains arriving af- 
ter Overlap- These chains will have / values which decrease 
with time since some monomers will be unable to bound 
due to the presence of earlier arrivers. This would lead to 
a continuously broadening spectrum of / values with in- 
creasing time. In practice this overlap time is often close to 
the saturation time (e.g. for N = 1000, Overlap = -25f^ t em ) 
so this broadening effect has little time to develop. 



4.2 Late Stages: Diffuse Outer Layer 

As the surface density approaches saturation, a /bound ~ 
1, the availability of surface sites on the surface becomes 
scarce and the late coming chains cannot fully adsorb on 
the surface as the early arrivers did. Let us suppose that 
the establishment of one surface attachment requires an 
empty spot large enough to accommodate n con t bound 
monomers, where n con t is a small-scale species-dependent 
number similarly to u). The surface density of these "super- 
sites," Psupor, is becoming smaller with time and the mean 
separation, Z scp , between neighboring supersites increases 

— 1/2 

accordingly, Z sep w p S uper- Now in order for late-coming 
chains to adsorb onto these surface spots, they have to 
form loops joining up these sites as shown in fig. [!|] We 
model the adsorption of chains at these late stages by as- 
suming that the size s of such loops is the equilibrium 
subcoil size corresponding to l sep , i. e. 



fl „3/5 _ , _ -1/2 
u-0 — l S ep — Fsupcr 



(36) 



Thus chains which adsorb at the instant when the typ- 
ical loop size in eq. I|36|l is s, have a fraction of bound 
monomers given by 



dr hc 



or 



(37) 



Now as /bound approaches its asymptotic value, rj°£ und , 
which is another nonuniversal small-scale dependent quan- 
tity of order a -2 , then 



An 



bound ^contPsuper — ^cont g^g 3 



(38) 
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Fig. 10. Sketch of predicted final layer structure resulting from 
irreversible polymer adsorption. The layer consists of two parts 
(one chain from each part is highlighted): (i) An inner region of 
flattened down chains making u>N contacts per chain, where to 
is of order unity, (ii) A diffuse outer layer build up from chains 
each making fN <C N contacts with the surface. The values of 
/ follow a distribution P(f) ~ / -4 ^ 5 . Each / value corresponds 
to a characteristic loop size for a given chain, s ~ n CO nt//. 

where ^r boU nd = ^bSund ~ Abound and we used eq. (jSEJ. 
Eq. Q38[) is simpler to understand starting from the com- 
pletely saturated surface: it states that if one were to un- 
peel chains from a completely saturated surface so as to 
create p SU pcr supersites, the number of bound monomers 
freed would be n cont per supersite. 

Now integrating eq. (|37|1 after using eq. Ij38|> one has 

a ^^bound = ^cont(o Ar/6) 6 , (39) 

which together with eq. (|35|l of the early stages describe 
the full theoretical relation between bound and total frac- 
tions plotted in fig. Eta). Thus eq. (|39J) predicts a very 
sharp saturation of the bound fraction as r — > r°°. 

Now since in our picture every point along the /bound (-0 
curve corresponds to a unique value of /, see eq. IpTfjl. then 
the weighting of / values is given by 

p (f) = jk;J dr6 (f- dr hound /dr) 

= ^(\d 2 r houad /dr 2 \ arhomd/ar ^y 1 , (40) 

which after using eq. i|39|) . leads to 

p(f) - cr A,b (f « i) (4i) 

withC«6/(5a 2 r oo [n cont ] 1 /5). 

The full theoretical prediction for the final layer's dis- 
tribution of / values is thus the sum of eq. I|41l) and a 
sharply peaked function at u> from the early stages. The 
overall distribution is shown in fig- Elb) , binned into bins 
of width Af = 0.02. The binned distribution exhibits 
2 peaks, of different origin: the peak at large / corre- 
sponds to early arriving chains and its position is species- 
dependent while the peak at / = represents a diffuse 
outer layer shown in fig. 1101 and is due to the universal 
small / form, P(f) ~ J" 4 / 5 . 

Let us consider now the distribution of loop sizes in the 
resulting diffuse layer, f2(s), equal to the number of loops 
of length s per unit area per unit loop length. This may be 
found by noting that the bound mass corresponding to a 
certain <irb oun d has a unique s value and must thus equal 
the bound mass in loops of length s, namely n con t^2(s)ds. 



Using eq. jSHjl this leads to 

f2(a) « a- 2 s- n ^ . (42) 

What density profile c(z) does this loop distribution gen- 
erate as a function of distance z from the surface? To 
determine the density profile we follow similar arguments 
to those of refs. [T^lfT^lfT^] which relate loop distributions 
to density profiles in both specific and general cases. We 
do that by assuming that the density at a given height z is 
due to contributions from loops longer than <j(z), where a 
is the loop length which extends spatially up to height z. 
The density at a given height is determined by the number 
of loops which are long enough to reach this height: 

c ( z ) « f ds f2( s ) . (43) 

Ja(z) 

Since in our case loops are not stretched one has z ~ 
acr 3 / 5 . Thus using eq. IpT^I in eq. (|4*3")l one obtains the 
algebraic profile, 

/OA 4 /3 

c(z) « c(a) (-) , (44) 

which interestingly has the same scaling form as de Gennes' 
self-similar profile of equilibrium polymer layers |10j . 

Let us finally consider the kinetics of the building up 
of the diffuse layer. Since the rate of attachment is pro- 
portional to the density of available surface sites, the early 
kinetics of eq. (|33|l generalize to 

dr 

-jr = 4>buikQAr bound . (45) 

(In eq. (|45|l we did not include free energy barriers due 
to loops of the already partially formed layer which would 
present exclude volume repulsion to new chains arriving 
at an empty supersite. In fact, in order for our picture to 
be self-consistent such effects must very small; were there 
any dangling loops near an empty supersite they would 
adsorb onto this site.) Thus from eq. (|45(l one has 

a 2 Ar^F(Cr/t) 1/5 , 

a 2 AT bound » G (Cr/t) e/5 , (t > tfm (46) 

after using eq. The prefactors F — 6(5n cont /6) -1 / 5 , 
G = «cont(C/6) 6 are close to unity. Together with the 
short time kinetics, eqs. (|33|l and H34[) ■ these evolutions 
are sketched in fig. ITTI 

So far in this section good solvent conditions were as- 
sumed. Generalizing to the case of theta solvents is straight- 
forward, by replacing v — 3/5 — > 1/2. One finds quali- 
tatively similar results with zi_Tb ounc i ~ e~ Fa replacing 
eq. I|39|l . and the distribution of bound fractions is now 
P(f) ~ replacing eq. (|4*T|l . 

5 Irreversible Physisorption 

We consider now the other important class of irreversible 
adsorption: strong physisorption. In this case monomer- 
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Fig. 11. Time evolutions of total, r, and surface bound, 
Abound, monomers per unit area, as a function of time, as 
predicted by theory. Both chemisorption (thin lines) and ph- 
ysisorption (thick lines) are sketched. The total mass for 
chemisorption grows initially linearly, r ~ t, and then slows 
down as the surface saturates, Ar = r°° — r ~ t" 1 ^ 5 . For 
physisorption F ~ t 1 / 2 has diffusion-controlled form for all 
times until saturation. The surface bound part for chemisorp- 
tion undergoes three regimes: initially /bound ~ t 8 ^ 3 , cross- 
ing over to /bound ~ t and then to the late stage satu- 
ration behavior, AFbornxd ~ t~ 6//5 . The surface bound part 
for physisorption follows initially diffusion-controlled kinetics, 
/bound ~ i 1 , but for longer times the kinetics slow down, 

AZ~bound = -/bound — /bound ~ [1 — {t/t^J/ ■ 

surface bonds form immediately upon contact and ad- 
sorption kinetics are diffusion-controlled. Despite the com- 
pletely different kinetics compared to chemisorption, we 
find the resulting polymer layers have nevertheless almost 
identical structure. 



5.1 Early Stages: Monolayer Formation 

Initially the surface is empty and it is thus inevitable that 
any chain whose center of gravity diffuses within the coil 
size [731I5T] R F = aN 3 / 5 of the surface will adsorb. Thus 
attachment of chains is diffusion-controlled and the sur- 
face coverage grows as 

a ^r(t) « ^ (Dt) 1 / 2 , (47) 
a 

where D is the chain center of gravity diffusivity. 

Consider now the part of adsorbed mass which has 
bonded with the surface, /bound- Immediately after the 
first attachment, the subsequent monomer arrivals on the 
surface occur at the rate of their diffusion on the surface. 
Imagine there were no reactions and that the surface was 
a penetrable plane. Then, given the coil is initially next to 
the surface, all monomers would cross the plane at least 



once within the bulk coil relaxation time Tb u ik |79II81| . 
With reactions turned on, each time a monomer reaches 
the plane, it would react with it. We expect the effect of 
the resulting constraint to further accelerate the rate of 
new monomer arrival onto the surface and thus the chain 
physisorption time r a d s would be at most T DU ik- We do not 
attempt to analyze the details of these kinetics involving 
polymer hydrodynamics with increasing number of con- 
straints. Such kinetics would anyway be very difficult to 
detect experimentally since typical coil relaxation times in 
solution are microseconds. Here, it is enough to know that 
Tads ~ Tbuik, which is supported by numerical simulations 
in the absence of hydrodynamic interactions in ref. |70| . 
Thus provided the solution is dilute, (f> < (ft*, each chain 
adsorbs fast enough onto the surface before other chains 
interfere with it. 

Since the single-chain adsorption time is faster than 
new chain arrival on the surface, a monolayer of flattened 
chains starts to develop on the surface, just as in the case 
of chemisorption. Defining w to be the surface-bound part 
of a completely collapsed chain, as in chemisorption, we 
thus have 

Abound = , (48) 

This is the initial linear regime in fig. ^a). 

The kinetics of eq. (|47() proceed up to a time at which 
the density of available surface sites starts to become small, 
i. e. when r rs a~ 2 . This occurs at a time of order 

C h t yS «r bulk (-f-) 2 . (49) 

V 0bulk / 

After this time, new chains must form loops joining dis- 
connected empty sites and a diffuse layer starts to build 
up. 

5.2 Late Stages: Diffuse Outer Layer 

For times longer than ig at ys , the late coming chains start 
to see a continuously decreasing density of available sites 
for adsorption. Will the adsorption kinetics continue to be 
diffusion-controlled? Consider a chain which was brought 
by diffusion to within one coil size Rp of the surface. Let 
us see if a bond forms before it diffuses away. The number 
of collisions the chain makes with a certain site on the wall 
within Rp of its center of gravity before it diffuses away 
is: 

Na 3 

N co \l S3 (TbulkAo) ^sitc , N sitc S3 r — — . (50) 

rip 

Here t a is the monomer relaxation time and AT S £ te is 
the mean number of the coil's monomers per surface site 
within the coil's projected surface area. This would sim- 
ply be of order Na 3 /Rp, but the presence of the hard wall 
reduces it by a factor r = Z sur f(N)/Zb u ik(N) = 1/N as 
in eq. Using [791I5T] Tbuik ~ t a (Rp/a) 3 one has from 
eq. J5T5J) 

JVcou w 1 (51) 
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Thus the chain makes of order one surface contact with 
every site on the surface within the coil size. It follows 
that unless the surface density of available sites is so small 
as to be of order one per R"p, it is inevitable that dur- 
ing the chain's residence time on the surface at least one 
monomer-surface bond will form. Thus the diffusion - con- 
trolled kinetics of eq. 1(47(1 continue up to the point where 
essentially every empty site is filled. 

What is the mode of adsorption of these late coming 
chains which attach to the few empty sites? We model 
adsorption in these late stages in a similar way to the 
build up of the diffuse outer layer in chemisorption (sec- 
tion 4.2). The late coming chains, after adsorbing on an 
empty site large enough to accommodate n con t monomers, 
form bridges to nearby empty sites which are loops of s 
monomers when the average separation between neighbor- 
ing empty sites is / sop « as" (see fig- EJ). This separation 
becomes larger as more and more chains adsorb. Thus the 
diffuse layer density profile c(z), distributions of loop sizes, 
/2(s), and distribution of fraction of adsorbed monomers 
per chain, P(f), are the same as for chemisorption (sec- 
tion 3) . Including the monolayer contribution due to short 
times, the resultant overall /bound (-0 and P(f) relations 
are plotted in fig. This is a repeat of the chemisorp- 
tion curves with the exception of the early /bound ~ r 8 ' 3 
regime exclusive to chemisorption. The sketch of the re- 
sulting polymer layer is the same as fig. ITU1 

The major distinction between chemisorption and ph- 
ysisorption arises in the kinetics of layer build up. The 
difference in the r(t) kinetics is evident by comparing eq. 
I|47l) for physisorption, valid up to complete surface satu- 
ration, to eqs. I|33[l and (|46|l for chemisorption. Similarly, 
the kinetics of the bound fraction, /bound are also very dif- 
ferent. For physisorption, the short-time kinetics are given 
by substituting eq. (|47|l in eq. 148(1 while the long-time be- 
havior is found using eq. (|47|l in eq. (|39|l . leading to 



(t « & 



^bulk» 1 {Dt) 1 / 2 ^ V^"finalJ, r) 



Here tg^i * s * ne complete surface saturation time, af- 
ter which the diffuse layer has completely formed. This 
timescale is of the same order of magnitude as t^l s since 
the surface coverage of the outer layer, J F dz c(z) « a~ 2 , 
is of the same order as the monolayer coverage. Thus from 
eq. (|47|l one sees that the time for diffuse layer formation 
is of the same order as the time of monolayer formation 
since both require the diffusion to the surface of about the 
same quantity of mass. The complete evolution of r(t) and 
Abound (t) are sketched in fig.HTl 



6 Discussion 

6.1 Comparison of Theory with Experiment 

We studied theoretically the structure of polymer layers 
formed by irreversible adsorption onto surfaces from di- 



lute solution under good or theta solvent conditions. Rel- 
atively few theoretical works have addressed irreversibil- 
ity in polymer adsorption. In refs. |67II68| irreversible ph- 
ysisorption was studied analytically and numerically within 
the framework of self-consistent mean field theory. These 
workers found that compared to equilibrium, the irreversibly 
formed layers are different in that (i) the asymptotic sur- 
face coverage r°° was larger in equilibrium, and (ii) in the 
irreversible case the density profile c(z) was found to be 
larger for small and large z but smaller for intermediate 
values as compared to equilibrium. In another work |4b) 
single chain chemisorption of PMMA onto Al was studied 
by solving numerically the loop kinetics, eq. (fP3|) . using 
chain statistics corresponding to theta-solvent solutions. 
It was found that the chain adsorbs in a zipping mecha- 
nism the origin of which was greatly enhanced reactivities 
for monomers neighboring a graft point. These reactivities 
were taken as greatly enhanced on the basis of electronic 
structure calculations [82U83) . 

In this work the cases of irreversible physisorption and 
chemisorption were examined separately. We found that 
in both cases the final layer consists of two regions. In 
an inner region of completely flattened chains each chain 
has on average uiN monomers bound to the surface. The 
value of lu is species-dependent and represents the effect 
of steric constraints at the monomer level which prevent 
all monomers of the chain from bonding. The outer re- 
gion is a tenuously attached diffuse layer of chains making 
fN <C N contacts with the surface. The distribution of / 
values among chains is universal and in a good solvent fol- 
lows P(f ) ~ ./"~ 4 / 5 for / <C 1. This double layer structure 
is illustrated in fig. 1101 which exhibits two peaks at small 
and large /, respectively. Our theory cannot capture the 
numerical value of the peak amplitude near f = u> but 
the existence of a peak (or possibly of a depletion region) 
is due to small-scale effects. The first layer of flattened 
chains is special because these arrive at a bare surface, 
and chain features on the scale of a monomer size are in- 
volved. In our picture, one can think of the next set of 
arriving chains (which are just starting to build up the 
outer tenuously attached region) as seeing a different sur- 
face, one containing a certain areal fraction of empty sites 
less than unity. This set cannot flatten down completely. 
The third set of arriving chains sees a surface containing 
a yet lower density of free sites, and so on. For large layer 
numbers, a chain sees a universal surface whose features 
evolve self-similarly with increasing layer number. This 
leads to a universal power law form for P(f). Returning 
to the earliest layers, these are formed before universality 
onsets, reflected in a non-universal feature in P(f) near 
f = u. 

Using infrared absorption spectroscopy, measurements 
of P(f) for irreversible physisorption of PMMA onto ox- 
idized silicon have been pioneered in the experiments of 
refs. |33H31ll3*2] measurements from which are reproduced 
in fig.l^b). In this figure, the theoretical prediction is com- 
pared to experiment by binning values in ranges Af — 
0.02 and converting analytical predictions to a histogram. 
Since / = 0.08 is the lowest observed experimental value, 
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we cut off the theoreticaf distribution at the same point 
for the sake of comparison. The contribution from the in- 
ner layer is represented by a delta-function centered at 
lu fti 0.47, though as discussed in section 4, we expect 
this to be broadened. We see that the agreement with ex- 
periment is good and captures the bimodal shape of the 
distribution. We note that measuring lower / values to see 
a clear signature of a power law regime may be difficult: 
our theory suggests these chains are a small fraction of the 
total chain population (the integral of P(f) is dominated 
by large /). However since those chains lie in the outer- 
most region of the layer, they would determine important 
physical properties such as hydrodynamic thickness and 
the strength of interaction of the polymer layer with an 
approaching interface. In comparing with theory above, 
we have used our predictions for good solvents. However, 
the experiments of refs. |31U32j were performed under sol- 
vent conditions slightly better than theta. Thus it may 
be that a more appropriate comparison is with our the- 
ory specialized to theta solvents, for which the predicted 
P(f) ~ 1// is very similar. 

The building up of the double layer is apparent in the 
shape of Ibound(-T) where /bound, r , are surface-bound 
and total surface coverage, respectively. For irreversible 
physisorption we found that initially /bound ~ wT. As 
the surface saturates, the diffuse layer starts to develop 
and ( r b°ound " Abound) ~ (^°° - rf as the asymptotic 
values (denoted by oo) are approached which are of or- 
der a~ 2 , where a is monomer size. We thus predict a very 
sharp saturation of the bound fraction as r saturates; 
the resulting curve is plotted in fig. Ufa). A very simi- 
lar curve has been measured in the experiments discussed 
in the previous paragraph |31II32| . where the initial slope 
was u w 0.47, the value we used in fig. |^b). In addi- 
tion the shape of the curve was found to be independent 
of chain length, which is also consistent with our model. 
For chemisorption, the curve is similar but we predict an 
additional /bound ~ r 8 / 3 regime as shown in fig. Efa). 

Contrary to physisorption, the experimental picture 
for chemisorption is much less clear, though since time- 
scales are intrinsically much longer, the kinetics might 
be easier to probe. Experiments on various systems have 

been performed |EI|13|EllEHllEnillE| the results of which 
however cannot be interpreted within the framework of 
our present theory since they involved simultaneous ph- 
ysisorption and chemisorption. In certain cases the degree 
of polymer functionalization was also varied HSUHHEU 
Clearly, the dependence of resulting polymer structures 
on degree of polymer or surface functionalization is an 
interesting aspect deserving further study. 

This paper has addressed both the structure of the 
final layer and the irreversible kinetics of its formation. 
Chemisorption kinetics are very slow and are chemically 
controlled. We found that initially JT ~ t, followed by r ~ 
t^ 1 / 5 as the surface saturates. We remark that for long 
enough times, adsorption onto a planar surface always be- 
comes diffusion-controlled after a Q-dependent timescale 
|S3ISB|i where Q is reaction rate upon monomer-surface 
contact. Here our assumption is that Q is small enough 



such that this cross-over time is very large, which is the 
typical case for ordinary chemical species. For physisorp- 
tion, diffusion controlled kinetics apply until saturation, 
r ~ t 1 / 2 . The corresponding time dependencies for /bound 
follow from /bound (r) an( i are shown in fig. 1111 

A large part of this work has focused on single-chain 
chemisorption. Unlike physisorption where single chain 
adsorption is complete at most after the coil's bulk re- 
laxation time, Tbuikj we found that the adsorption time in 
chemisorption is much larger, r a d s ~ Q" 1 ^/ 5 . Thus for 
typical Q w Is -1 , A" = 1000, one has r a d s « 60s, a time 
accessible to experimental measurements. During r a d s the 
chain adsorbs onto the surface in a mode we call accel- 
erated zipping: initially the chain adsorbs in a zipping 
mechanism growing outwards from the first attachment 
point with the surface, but with increasing time distant 
monomers adsorb forming large loops and new sources 
for further zipping which accelerate the chain's collapse. 
The distribution of these loop sizes follows a power law 
w(s) ~ s -7 / 5 while the number of bound monomers grows 
as 7bound(i) ~ t 5 / 3 . 



6.2 Differences between Irreversible and Equilibrium 
Layers 

We conclude with a general comparison between the fi- 
nal non-equilibrium layers predicted by our theory and 
equilibrium layers. This discussion is limited to good sol- 
vents. Our results for the final loop distribution of the 
layer, i?(s) ~ s~ n / 5 and density profile c{z) ~ z~ 4 / 3 
have interestingly identical scaling form to the correspond- 
ing equilibrium results. On the other hand, the configura- 
tions of individual chains are very different. In the equi- 
librium layer a given chain has ND(s) loops of length 
s or greater, where D(s) = ds' Q{s') ~ s~ 6 / 5 . Be- 
cause of screening effects these are essentially independent 
blobs and their 2D spatial extent parallel to the surface 
is [ND(s)] 1 / 2 as 3 / 5 = aN 1 / 2 . This is true for all scales 
s; in particular, there is order one loop of length AT 5 / 6 , 
also of size aN 1 / 2 . Hence a typical chain has a size 
of order aN 1 / 2 , the ideal result (to within logarithmic 
corrections Chains in the non-equilibrium layer are 
by contrast of size aN 3 / 5 , due to the adsorption kinetics 
which in the case of chemisorption entailed occasional ad- 
sorption of very large loops of length A" which ensured the 
final state of the chain extends a distance of order aN 3 / 5 
parallel to the surface. We did not attempt to unravel the 
complex details of physisorption, but our assumption was 
that the 3D coil ends up being roughly projected onto the 
2D surface, leading again to a lateral size aN 3 / 5 . 

The most fundamental distinguishing feature of the 
non-equilibrium layer is that different chains have differ- 
ent loop distributions and fractions of adsorbed segments, 
/, which are frozen in time. There is no longer one sin- 
gle thermodynamic /. Instead there are infinitely many 
classes of chains, each with its own / value, and the num- 
ber of chains in each class is proportional to P(f). Each 
chain has a characteristic loop size, s ~ 1//. 



Ben O'Shaughnessy, Dimitrios Vavylonis: Irreversible Adsorption from Dilute Polymer Solutions 



15 



This is very different to the equilibrium layer where 
all chains are statistically identical, in that every chain 
explores in time the same range of / values with the same 
relative weighting P cq (f). In fact, one can show that for 
large N this distribution P cq is very sharply peaked at a 
mean value / of order unity, with width of order TV -1 / 6 . 
(The mean value is due to small loops of order unity while 
the fluctuations are due to the mass in very large loops 
with length of order iV 5 / 6 . Since there is order one such 
loop, this diminishes / by an amount of order iV~ 1//6 .) 
Thus even at any fixed moment in time it is still a true 
statement that almost all chains have the same / value, 
/, to within fluctuations which are small for large enough 
N. We remark that there is a small population of chains, 
a fraction of order N~ x / 5 of the total, with / values far 
removed from the mean: these are the chains with tails |15j 
of length of order N which determine the layer height 
w aiV 3 / 5 . 

These differences in chain statistics between irreversible 
and equilibrium layers have important implications for 
various physical properties of the layer. Consider a ph- 
ysisorbing system whose surface relaxation kinetics are 
very slow, but not truly irreversible (we expect this is a 
rather common case). In this case for times of order the 
layer formation time we expect the layer is well-described 
by the irreversible P{f). For much longer timcscalcs on 
which surface bonds are reversible equilibration processes 
will slowly evolve P(f) into P eq (f). Our picture predicts 
that during this process the total coverage, density profile, 
and loop distribution remain unchanged to within pref- 
actors of order unity. Single chain statistics will be sig- 
nificantly affected, however, and individual chains must 
shrink in order to occupy a much smaller region in space. 
This change is expected to profoundly modify the physical 
properties of the layer. For example, the fraction of chains 
with long loops and tails of order N will be reduced by a 
factor N^ 1 / 5 , thus strengthening the outer region of the 
layer which is exposed to the bulk. Similarly in experi- 
ments probing the kinetics of exchange between chains in 
the bulk and those in the layer |21l251l351l27ll25lli&Ha01 
I3lfll3l)ll3"7ll38ll3"§] . we expect strong aging effects since the 
number of easily exchangeable, loosely attached chains is 
decreasing with increasing aging time. This would lead to 
a decrease of the initial exchange rate with aging time as 
observed 28 , 29 37 38 . Further theoretical work is clearly 
needed to interpret the phenomenology of these exper- 
iments. One unexplained observation, for example, con- 
cerns poly(ethylene oxide) adsorption onto silica where 
the bulk-surface exchange rate of chain subpopulations in 
non-equilibrium layers has been measured to be indepen- 
dent of the time the subpopulation was incorporated in 
the layer 

What is the origin of the theoretically predicted irre- 
versible layer's overall loop distribution being the same 
as that in equilibrium? This is rooted in the kinetics of 
our model. In the late stages, we assumed the equilibrium 
relation between loop size and number of monomers s in 
the loop, R ~ s v . That is, a given loop is assumed to 
have equilibrium statistics. It is important to note that 



this does not by itself lead to the equilibrium loop distri- 
bution. Our basic assumption on the late stage kinetics is 
that the kinetically selected loop size R at any moment 
is determined by the current density of free supersites, 
Pauper ~ ~L/R 2 - Taken together with R ~ s" , this then 
leads to the S2(s) of eq. I|42(l. These kinetics are mean 
field in character, assuming uniformly smeared free sites 
determined by the current global value of p supcr rather 
than any local feature. The correct choice of kinetics is 
governed by the kinetics of the chain adsorption process. 

To illustrate this point, consider a general situation, 
involving chain adsorption onto a cZ-dimensional surface. 
Suppose the mechanism of adsorption is pure zipping and 
we follow the process from a starting situation where the 
surface has a uniform density of empty surface sites, /9 S upcr : 
l~ d . Roughly speaking such a zipping chain performs a 
simple random walk on the surface, each step filling one 
empty site and producing a loop of length s = l 1 /" '. The 
statistics of this process depend on dimensionality, (i) Sup- 
pose d = 1. Since a ID random walker explores space 
"compactly" [HH], i-e. visits each site within its exploration 
volume many times, by the time the zipping is complete 
and the chain has fully adsorbed a surface patch depleted 
of empty sites will have been created in the region where 
the chain has adsorbed. Later-arriving chains can then 
occupy only the empty regions between such patches cre- 
ated by previously adsorbed chains. We expect the den- 
sity of free sites to fluctuate very strongly, invalidating any 
mean field approach, (ii) d=3. The zipping process is now 
a "non-compact" random walk, occupying only a small 
fraction of the free surface sites within its exploration vol- 
ume. Fluctuations are small, and the mean field approach 
which assumes an essentially uniform distribution of free 
sites of density p S upcr decreasing quasistatically in time is 
valid. This is the assumption implicit in the present pa- 
per, (iii) d=2. Random walks are now marginally compact 
and one expects logarithmic corrections to the mean field 
picture. 

The situation treated in this paper is of course d = 2. 
Thus one might expect logarithmic corrections to the loop 
distribution of eq. (|42|l which was arrived at via the mean 
field argument. In fact, for two reasons we expect the 
zipping process to execute a surface walk somewhat ex- 
panded from a simple marginally compact 2D random 
walk. Firstly, no site can be adsorbed onto twice by the 
zipping random walk. Thus walks are repelled from fre- 
quently visited areas. This is a kind of self-avoidance, and 
one anticipates swelling of the walk somewhat similarly to 
the "true" self-avoiding walk jHD] and the "kinetic growth" 
walk . Secondly, the chain adsorption processes we 
have studied are not simple zipping. For chemisorption, 
the zipping is accelerated by occasional very large steps 
producing loops of all sizes up to N. These large excur- 
sions result in the flattened chain occupying an area of 
order N 6 / 5 rather than N. Since the fraction 7V~ 1//5 of 
occupied sites within this area is small, one returns to a 
non-compact situation where mean field applies. In reality, 
this is over-simplified since local zipping from nucleation 
points presumably involves correlations of the type charac- 
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terizing pure zipping, so non-uniform depletion of empty 
sites may occur near zipping centers. For physisorption 
the process is even more complex, presumably involving a 
combination of local diffusion- mediated zipping and large 
loop formation, but one expects a projection area onto the 
surface again of order TV 6 / 5 which naively suggests non- 
compact character. 

From the above remarks, it is clear that stepping be- 
yond our rather simple approach to irreversible adsorption 
involves many complex issues. We hope our work moti- 
vates future experiment and theory in this direction. 
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under grant no. DMR-9816374. 



and to x = s' / s in the second term, on the rhs of eq. 1)19(1 . 
respectively, one has 

ds' w T {s')k(s\s') — I ds' to T (s) k(s'\s) = 



2 f 1 , ss\ CO) u>t(s) f 1 



dxLU T ( —] — : 

s" ,/ \xJ x 1 -" 



dx C(x) (54) 



after using eq. I(ll|l . Multiplying eq. ((54(1 by s, integrat- 
ing over all positive s and changing variables from s to 
y = s/x, one finds that the result vanishes by virtue of 
the identity 2 J Q dxx£(x) — J Q dx((x) for ((x) symmet- 
ric around x = 1/2. Thus M tot is conserved. 



A Single Chain Chemisorption: Interior 
Monomers Attach First 

In this Appendix we show that a single chain is much 
more likely to chemisorb on a surface by one of its inte- 
rior monomers rather than with one of its ends. Consider 
the s th monomer, where by definition s is the length of 
shorter part of the chain, the other having length N — s. 
The reaction rate of this monomer with the surface is pro- 
portional to the partition function of the chain anchored 
by this monomer on the surface: 



■^cont(s) 
7]{x) 



li N r]{s/N)N 



72 — 1 



1 



(x < 1) 
(x - 1/2) 



(53) 



where fi is a constant of order unity. We derived Z con t(s) 
by demanding that (i) it is a power in s and (ii) in the lim- 
its s — > 1 and s — > N/2, one obtains the known partition 
functions Z cont (l) « fi N N~ <1 ~ 1 , Z cont (N/2) a ^ N N^~\ 
respectively. The numerical values of 71 , 72 are |77ll?)ll7Hj 
71 ~ 0.68, 72 =7—1 w 0.16, where 7 is the susceptibility 
exponent [75| . 

Now the total rate for any monomer to react is pro- 
portional to J^ 2 dsZ cont (s). Since 72—71 > —1 this 
integral is dominated by large s. Hence it is much more 
likely for a monomer in the middle to be the first one to 
react with the surface, rather than an end, roughly by 
a factor NZ cont (N)/Z cont (l) « N 1+ ^~^ « A^ 48 . Thus 
the accelerated zipping always propagates outwards from 
an interior monomer, as schematically illustrated in fig. El 



B Proof that Eq. (JTH|) Satisfies Mass 
Conservation 

We show in this Appendix that eq. 1(19(1 conserves the total 
number of monomers M tot = J Q ds s w T (s), by showing 
that dM tot /dr — 0. This is equivalent to showing that the 
integral of the rhs of eq. ((191) over all positive s is zero. 
Changing variables from s' to x = s/s' in the first term, 



C Derivation of Eq. ((231) from E q- 

Defining £ = N — s and w T (£) = u T {N — £), one has from 
eq. (ESI) 



(55) 



Laplace transforming £ — > E (we allow £ to take any 
positive value, i. e. we seek the A^ — * 00 solution of eq. 
((23(1 ) one has from eq. I(55|) : 



— L) T {E) = -AQE v Zj t {E) , 
dr 

which after integration over r becomes 
u T (E) = e- A Q TE " . 



(56) 



(57) 



We used the initial condition uio{E) = 1, since wo(£) = 
5(£) (see eq. (J23J)). Laplace inverting eq- lf57|) one re- 
covers eq. (J23J. 



D Steady State Solution of Eq. (I?T|> 

We seek a quasi-static power law solution, uj t (s) = H/s a , 
to eq. l(2Tll . Substituting to the left hand side of eq. (|77jl 
one has 



ds' \{s\s')uj t (s') = 
QH 



1 r l-y 
v-2+a 



dyy" 



dx((x) , (58) 



after changing variables to y — s/s'. Here we took the 
well-defined limit A^ — > 00 since the solution we seek 
should be independent of the tail's length. The only value 
of a for which eq. I(58|) is zero is a = 2 — v. The self- 
consistency of the quasi-static solution, eq. 1(28(1 . can be 
verified by substitution in eq. 1(21(1 and showing that for 
loop sizes where the solution is valid (s < s'° ax ) correc- 
tions arising from integration over s' > s" lax are small. 
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